Entropically driven transition to a liquid-crystalline polymer globule 
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A self-consistent-field theory (SCFT) in the grand canonical ensemble formulation is used to 
study transitions in a helix-coil multiblock copolymer globule. The helices are modeled as stiff 
rods. In addition to the established coil-globule transition we show for the first time that, even 
without explicit rod-rod alignment interaction, the system undergoes a transition to a nematic 
liquid-crystalline (LC) globular state. The LC-globule formation is driven by the hydrophobic 
helical segment attraction and the anisotropy of the globule surface energy. The full phase diagram 
of the copolymer was calculated. It discriminates between an open chain, amorphous globule and 
LC-globule. This model provides a relatively simple example of the interplay between secondary and 
tertiary structures in homopolypeptides. Moreover, it gives a simple explanation for the formation 
of helix bundles in certain globular proteins. 

PACS numbers: 61.30.Vx (Polymer liquid crystals), 87.14.Ee (Proteins), 87.15.-v (Biomolecules: structure 
and physical properties) 



The formation of secondary structure in proteins and 
in particular the a-helix-coil transition in homopolypep- 
tide chains is one of the well-investigated conformational 
transitions in biomolecules. In the well-known Zimm- 
Bragg (ZB)-theory of this phenomenon the polypep- 
tide molecule is considered as a one-dimensional coopera- 
tive system and the problem can be solved exactly. How- 
ever, ZB-theory is only valid in a fully denaturated state 
when all three dimensional interactions are negligible. In 
real systems the helical parts are often hydrophobic, such 
that this hydrophobicity drives the helix-coil copolymer 
into a globular phase. Furthermore, helical parts can be 
seen as rigid rods and hence a rod-rod alignment inter- 
action could be taken into account. Considering all these 
facts, it is natural to pose the question: how does the 
chain compaction affect the secondary and tertiary struc- 
ture? This problem has been partially discussed within 
computer simulations of globular proteins 0, . Some 
preliminary theoretical results concerning the formation 
of a LC-globule can be found in the review For a 
system of fixed composition of stiff parts (helices) and 
flexible parts, micelle formation has been studied using 
scaling considerations |^. 

In this paper the phase behavior of a helix-coil copoly- 
mer is discussed by means of a self-consistent field theory 
(SCFT). In addition to two- and three-body contact in- 
teractions between the segments, each segment is also 



allowed to undergo a microscopic transition from flex- 
ible to stiff (i.e. being part of a helix). To describe 
this local helix-coil transition a grand canonical formal- 
ism was used. A formal derivation of the grand canoni- 
cal SCFT free energy functional F[ip{r), tp{r, u); yit, e, a] is 
given in ||7|. We introduced a field (p{r) associated with 
a fiexible segment at position r and a field ip{r, u) as- 
sociated with a helical (rod-like) segment at position r 
with orientation u. The local densities of fiexible and 
helical segments are given by Pc(r) = l/2ip{r)'^ and 
Ph(r) = 1/2 / (i^w.'0(r, u)^ respectively. It is very difficult 
to deal with the full orientation dependence of tjj{r,u.). 
However, the helices are modeled as stiff rods without chi- 
rality, so we expect the solution for ip{r) and tp{r, u) to 
have an uniaxial symmetry. Therefore we choose, with- 
out loss of generality, the z-axis as the preferred orienta- 
tion and expand the orientation dependence in terms of 
Legendre polynomials. To lowest non-trivial order this 
expansion reads 
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Using this expansion u can be integrated out and the 
SCFT free energy functional of our model is then given 
by the following form 
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The numerical prefactors in Eq. (O are due to the expan- 
sion of the orientation dependence - see Eq. - and the 
subsequent integration over u. The chemical potential fi 
is used to control the total number of segments of the 
polymer N. To model the energy gain due to the for- 
mation of hydrogen bonds in helices we introduced an 
energy gain per helical segment e. The cooperativity ef- 
fect in the formation of helices is taken into account by 
the cooperativity parameter a. It can be regarded as a 
fugacity of the interfaces (or junction points) between a 
helix and a flexible part, a — I means no cooperativity 
effect and a — reflects total cooperativity (i.e. the sys- 
tem can only form a fully flexible chain or one long helix) . 
The interaction parameters v and w are global two- and 
three-body interaction constants between all segments, 
whilst X controls the strength of a selective two-body 
interaction (due to hydrophobicity) between the helical 
segments only. Further technical details can be found in 


At this point it is important to emphasize that all inter- 
actions in this model are point contact interactions. We 
do not take into account any angle dependent interac- 
tions between the helical segments which explicitly favor 
alignment of two helices. The total number of segments 
N is given by the sum of all segments in the helical and 
coil state: N — A^h + Nc, where 

^ II d\ [V.o'(r)-HV'|(r)] 

Nc ^ \l dW{r). (3) 

Minimization of the functional in Eq. (0) with respect to 
if, ipo and yields a set of three coupled differential 
equations for the three self-consistent fields. These dif- 
ferential equations are highly non-linear and can only be 
solved numerically. This was done using the finite ele- 
ment toolkit Gascoigne 0. 

The system shows a collapse transition from an open 
chain to a dense globule. In this letter we always con- 
sider finite systems, therefore this transition is a crossover 
transition with a finite broadness. Hence we have to de- 
fine a point during the crossover as the transition point. 
Before we explain how we do this for the coil-globule 
transition it is necessary to clarify how we generally deal 
with the chemical potential. In the grand canonical en- 
semble the number of particles, here the total number 
of segments of the polymer N, is not fixed but its mean 
value is determined by equilibrium conditions. However, 
in a real experiment the helix-coil copolymer has a fixed 
length. In order to ensure this fixed length N we tune the 
chemical potential n for each set of physical parameters 
(w, w, X, e, cr) such that the equilibrium value of N, cal- 
culated by performing the integrations in Eq. Q is equal 
to the desired one. For a given set of parameters iV(/i) 
can be computed and a typical example of this is shown 
in FigH For /i ^ the total number of segments N 
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FIG. 1: N as function of for «; = 1.0, v = 0.2 and a = 10'*. 
The dotted curve corresponds to e = 0.08 and x ~ 0. The 
continuous curve corresponds to e = 0.1 and x = 0.0138. 



diverges. This corresponds to the N ^ behavior of a 
0-solvent chain. The divergence of TV at a specific value 
fi on the right hand side of the minimum corresponds to 
a fully collapsed infinite globule. The minimum of N{v) 
is naturally associated with the transition point. We al- 
ways fix N for a given set of parameters by tuning ^, so 
from a plot like the one shown in Fig^ we can distin- 
guish whether the system is left of the transition point 
(i.e. in the open chain regime) or right of the transi- 
tion point (i.e. in the globular regime). The transition 
point is now defined as the set of parameters at which 
the chosen fixed N is equal to the minimum of N{^). 
Both curves in Fig.(^ have their minimum at iV = 10^, 
the value which we have chosen for the results presented 
below. The minimum of the dotted curve corresponds 
to the first triangle (e = 0.08, x = 0) in the phase dia- 
gram Fig.® and the minimum of the continuous curve 
to the second triangle (e = 0.1, x = 0.0138). It is perti- 
nent to note that the characteristic dependence of N on 
fi for a simple homopolymer globule has first been dis- 
cussed by Kholodenko and Freed 0]. For all the results 
presented below, the global three-body interaction con- 
stant is set to w — 1.0, the global two-body interaction 
constant to v = —0.2 and the cooperativity parameter 
to cr = 10~*. The SOFT formalism is a reasonable ap- 
proximation only for systems in a globular state or close 
to a globular state. Therefore v has to be negative in 
order to ensure that this is always the case. Although we 
do not use the SCFT treatment of our model to describe 
an open chain of the helix-coil copolymer, it is perfectly 
possible with the definition given above to calculate the 
transition line between the globular state and the open 
chain state. 

The really astonishing feature of this model is that it 
also shows a crossover transition from a disordered amor- 
phous globule with low or mid fraction of helical segments 
©h = Nh/N to an ordered liquid-crystalline globule with 
very high fraction of helical segments. This transition 
occurs despite the fact that there is no explicit angle 
dependent alignment interactions. The transition is trig- 
gered by a subtle interplay of the entropy contribution 
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FIG. 2; Fraction of helical segments 0^ as a function of 
strength — x of the attractive two-body interaction between 
the helical segments only for different values of the energy 
gain per helical segment e. 
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FIG. 3: Nematic order parameter S as a function of — x for 
different values of e. 



(surface energy), represented by the derivative terms in 
Eq.©, and bulk interaction energy, represented by the 
X-term. It is well known jlQ] that in a simple homopoly- 
mer globule the surface energy has an entropic nature 
(since the conformational set of surface segments is con- 
strained) and is isotropic. As one can see from Eq. (0) hi 
our case the surface energy is anisotropic, so that after 
a proper inspection one can ensure that the surface 
tension in the xy-direction is smaller than the one in 
the z-direction. That is why the system tries to maxi- 
mize its lateral surface in xy-directions and minimize it 
in z-direction, i.e. a nematic, cigar shaped, LC-globule 
occurs. 

To measure orientational order in the system we define 
the nematic order parameter S as follows (see e.g. 0) 



J d^uP2{cos0)^/j^{r, u) 



1 

3iV 

^ / d^"^ Mr) [Mr) + VbMr) 



(4) 



Fig. ^ shows the increase in fraction of helical segments 
0h with Ixl during the transition and Fig.© shows the 
simultaneous onset of a finite nematic order parameter 
S. The onset of the transition is shifted to lower values 
of Ixl with increasing energy gain per helical segments 
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FIG. 4: The density of the flexible segments is shown on 
the left and the density of the helical segments on the right, 
e = 0.1 for all plots, x = -0.0138 in the top line, x = -0.0812 
in the middle line and x ~ —0.18 in the bottom line. 



e. This is due to an increase of bulk interaction energy 
for fixed x with increasing number of helical segments. 
The curves for e = 0.2 and e = 0.15 in Figs.|j2El start 
at non-zero values of x- These values of x correspond to 
the transition point between open chain state and glob- 
ular state, i.e. the triangles in Fig.©. The transition 
point for the crossover transition from an amorphous to 
a nematic LC-giobule has to be defined in a reasonable 
way. We define these transition points as the inficction 
points of the S'(x)-curves in Fig.©. 

To demonstrate how the shape of the globule changes 
during the transition we show in Fig.© a color-coded 
plot of the local density in g-z space, where g — 
y^x^ + y^. The center of the globule is in the bottom 
left corner, g is increasing from left to right and z is 
increasing from bottom to top. Red means high density 
and dark blue zero density. The top two pictures show 
the density profile at the transition point between open 
chain and amorphous globule. At this point the system is 
spherical and has a very broad surface layer of decaying 
density resembling the characteristics of a random walk. 
Although the density of the helical segments shown on 
the right is higher than the density of the flexible seg- 
ments shown on the left, their distribution and the shape 
of the profile is very similar. The middle two pictures 
show the density profile at the transition point between 
amorphous globule and LC-globule. The system adopts a 
slightly cylindrical shape indicating the onset of nematic 
order. It can also be seen that the density maximum of 
the flexible segments is not in the center of the globule 
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FIG. 5: Phase diagram of a helix-coil copolymer. The upper 
left area corresponds to an open chain, the lower left area to 
an amorphous globule and the right area to a nematic LC- 
globule. The little arrow to the left indicates the point in the 
phase diagram which corresponds to the top two pictures in 
Fig.©. The arrow in the middle corresponds to the middle 
two pictures and the arrow to the right to the bottom ones. 

denoting a repulsion of flexible segments from the center 
to the surface layer. The surface layer is now much nar- 
rower. The bottom two pictures show the density profile 
deep in the nematic LC-globule phase. The globule has 
developed a strongly asymmetric cylindrical shape indi- 
cating strong nematic order. The repulsion of flexible 
segments from the center towards the surface layer can 
be seen clearly and the surface layer is now very narrow. 

The results shown in Fig.Q allow us to compute a 
complete phase diagram of a helix-coil copolymer in e — % 
space, see Fig.©. The triangles are the transition points 
between open chain and globule. The squares are the 
transition points between amorphous globule and LC- 
globule. Note, the points plotted in the phase diagram, 
Fig.lO, are what we defined above as the transition 
points of rather broad crossover transitions. Therefore 
the boundaries in the phase diagram have to be under- 
stood as "center lines" of broader regions in which the 
crossover from one phase to the other occurs. Although 
the qualitative shape of the phase diagram stays the 
same, the position of the transition lines changes with 
A'^, a and v. If is increased the open chain region 
of the phase diagram becomes smaller (and disappears 
eventually) . With decreasing tr the transition from a dis- 
ordered to an ordered globule becomes sharper and the 
transition line between these two regions is shifted to 
smaller values of Since this transition occurs due to 
an interplay between surface energy and bulk interaction 
energy, the transition also becomes sharper for decreasing 
system size N , which, at first sight, is a rather unusual 
and surprising behavior. For smaller systems the surface 
energy plays a bigger role and therefore leads to a sharper 
transition. For N ^ oo the ordered globule phase finally 
disappears, since the surface contributions to the free en- 
ergy vanish for infinite systems. A detailed discussion of 



the modification of the phase diagram will be given in a 
subsequent publication [llj 

In summary, we presented a SCFT for a multiblock 
helix-coil copolymer chain based on grand canonical en- 
semble considerations. The minimization of the corre- 
sponding free energy functional - Eq. ^ - has been done 
numerically and it has been shown that three phase states 
can be clearly seen: open helix-coil chain, amorphous 
globule and nematic LC-globule. It is a novel result 
that the formation of a LC-globule occurs without ex- 
plicit alignment interactions between the helical parts. It 
is the entropical surface tension anisotropy which drives 
the globule in the nematic LC-state in order to maximize 
the density of the hydrophobic helical segments. In the 
presence of an explicit attractive alignment interaction 
(which can be taken into account by a Maier-Saupe term) 
one only sees an enhancement of this effect . We be- 
lieve that this model provides a relatively simple example 
of the interplay between secondary and tertiary structure 
in homopolypeptides. It can also give a simple explana- 
tion for the formation of helix bundles in certain glob- 
ular proteins. Both, simulations [l^ and experiments 
[l^ show that proteins can adopt not only the native 
state and completely denaturated state (open chain) but 
also so-called premolten and molten globular states. For 
helix-bundle proteins the premolten globule, which does 
not show any order of the helices, corresponds to our 
amorphous globule. The molten globule with ordered 
helices but without native contacts corresponds to our 
LC-globule. The chain length of proteins is typically of 
the order of A^ - 10^ - 10^ instead of A^ ~ 10^. However, 
as mentioned above the crossover from an amorphous 
globule to a nematic LC-globule with aligned helices be- 
comes even sharper for shorter chains This indicates 
that our model indeed provides a simple explanation for 
the formation of helix bundles. 
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